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Abstract:This paper presents the results of a study to develop a fast calculation routine which enables the 
prediction of the frost layer growth and the deposition velocity in terms of the surface temperature of the cold 
body, the rate of flow of the humid air and the rate of flow heat transfer and mass transfer equation for the flow 
of wet air around a circular cylinder are formulated and solved and the results are presented in function of the 
surface temperature, mass flow rate of the humid air and the mass flow rate of the working fluid. In order to 
validate the model and the numerical predictions, an experimental rig is constructed and adequately 
instrumented. Experiments were realized and the results were composed with the numerical predictions to 
confirm the validity of the model. 
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I. INTRODUCTION AND BIBLIOGRAPHIC REVIEW 

When humid air passes over cold surfaces the liquid particles are deposited over the surface in the form 
of crystals by complicated mechanisms of heat and mass transfer forming a porous layer. If the process of frost 
deposition continues for long periods this may lead to the formation of thick glacial layer of frost which impairs 
the performance of the equipment and increases the pressure losses. In these cases the formed frost may be 
removed by a variety of methods causing additional heat gains and the usual undesirable interruption of the 
cooling process. For these and other reasons it is necessary to study deeply the process of frost formation to be 
able to operate the equipments efficiently and for long periods. 

One of the first experimental studies on frost formation was realized by [1] in which the effects on the heat 
transfer and pressure drop were evaluated [2] estimated the effects of the frost growth on the energy losses in air 
conditioning and refrigeration systems and found that about 20% of the total energy is spent in the defrost 
process. Other important experimental contributions are due to [3], [4], [5], [6], [7] and [8]. 
Mau et al. [9] reported data and empirical correlations of the frost growth and its properties in terms of the 
temperature of the flat plate, air temperature, humidity, and velocity. 

Lee and Ro [10] studied the frost formation on the surface of a horizontal cylinder in a wet air stream and 
obtained results of the heat transfer flow rate around the cylinder, the effective thermal conductivity of the frost 
and the temperature distribution around the cylinder. Yang and Lee [11] proposed correlations for the frost 
deposition on the surface of the cylinder of the frost thickness, density, temperature of the surface of the frost 
layer. 

The literature survey shows that during the last fifty years many models were proposed to describe the 
phenomenon of frost formation as the model due to [3] which is strongly based upon experimental data and 
empirical equations. Parish and Sepsy [12] proposed a model to formulate the frost formation around a cylinder 
which by the use of global balances in the equations of momentum, energy and concentration try to quantify and 
evaluate the process of frost formation. The model proposed by Jones and Parker [13] uses global balances and 
considers the variation of the internal properties of the frost layer which makes the balances more representative. 
Hayashi et al.[14] elaborated a model which divides the process of frost formation into three successive periods, 
where the first is the one dimensional growth of the crystals, the second period is for three dimensional growth 
and the last period is the period of total development or nearly static growth. 

Tokura et al. [7] elaborated a set of parameters which enables the determination of the effective thermal 
conductivity. Aoki et al. [15] studied the third period of Hayashi et al. model in great details. 
Sami and Duong [16] proposed modifications in Jones and Parker model by using local frost properties to obtain 
better predictions. Padki et al. [17] proposed an approximate model for frost formation by the use of empirical 
correlations to estimate the heat and mass film coefficients. Sherif [18] presented a transient semi-empirical 
model for frost formation over a flat plate. The model estimates the surface temperature and the frost layer 
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growth using existent correlations for the heat transfer coefficient and Lewis' analogy to calculate the mass 
transfer coefficient. 

Tao et al. [19] proposed a frost formation model on a flat plate. The model is based upon the basic equations of 
heat and mass transfer and empirical correlations. The model solves the problem of frost formation in two 
stages, the first stage of short duration one dimensional crystal formation followed by the second stage which 
uses the data from the first stage as entry values. 

Scalon [20] solved the problem of frost formation over a vertical cylinder in a wet air stream. He used the 
equations developed by [19], but differently from Tao et al, he solved the flow field equations to determine the 
convection and diffusion coefficients between the flow and the surface of the frost layer. 

Salinas [21] and Ismail and Salinas [22] studied the case of frost formation over a flat plate placed parallel to the 
direction of wet air flow by using a formulation based upon the average volume local technique and were able to 
calculate the local density, frost layer thickness, mass deposition rate and phases mass fraction. They compared 
their results with available experimental and numerical data. Later Ismail and Salinas [22] extended their work 
to treat the case of flow of humid air and frost deposition between two parallel plates. 

Some recent experimental work was realized by [23] on the formation of frost on vertical cold cylinders of 
different arrangements including the case of two cylinders in series and three cylinders in a triangular 
geometrical arrangement. 

This paper presents the results of an investigation to develop a calculation routine which enables the 
prediction of the frost layer growth rate and the thickness of the frost layer under different operational 
conditions of mass flow rate of the working fluid, the surface temperature of the cold surface and the mass flow 
rate of the humid air. Experiments were realized on an experimental rig designed and instrumented to enable the 
variation of the air flow rate, the secondary working fluid flow rate and its temperature. The frost interface 
position was determined by a digital camera. The photographs were discretized to determine the real frost 
interface position. The experimental frost thickness and the frost deposition velocity were used to validate the 
numerical model and the numerical predictions. 

II. FORMULATION OF THE PROBLEM 

The mathematical model developed by [20] based upon [19] with some adaptation, was used in the 
present work to describe the process of frost formation around a vertical cold cylinder placed in a wet air stream. 
Accepting the model due to [19] ,at the moment is the most accepted model even it is not exact one, that the 
process of frost formation can be divided into two stages, where the first stage is one dimensional crystal growth 
followed by the second stage of interaction between neighboring crystals forming a porous matrix. Both Tao et 
al. and [20] used an empirical parameter called the transition time taken as 2/3 of the transition time given by 
[14]. The transition time is a parameter which indicates the transition predominance from one phase to the other, 
given in Fig. 1. 




258 280 2S2 264 266 

Temperature [K] 



Fig.l Values of transition times in the range of -20°C the 0°C. 
1.1 One dimensional crystal growth 

In this stage ice nucleus are formed on the cold surface and grow in the vertical direction as a consequence of 
the heat and mass transfer. 
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Following the model proposed by [19], the crystal growth is normal to the surface as shows in Fig. 2. Applying 
the energy equation on the element shown in Fig. 2. 
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Fig. 2 Physical model for the growth of the ice crystals, [19] 
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(i) 



(2) 



(3) 



Applying the mass diffusion equation one can obtain 



Pa 7Tt( d + d i dt f ~ d2 } dr = h m(nddr)(W r - W p ) 



p g -2d d ^-dr = h m (nddr)(W Y - W p ) 



(4) 
(5) 



Where Wp indicates the concentration evaluated at the temperature Tp of the solid phase and the result is: 
p g % = 2h m (W r -Wp) (6) 

These two equations represent the model for the first stage. The initial and boundary conditions necessary to 
solve this model are 

When r = 8 



d = d (] 



66 



dS h„ 



dt 



Pg 



(7) 
(8) 

(9) 
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When r = 0 



^ = 0 

dt 



When t = 0 
d = d Q 
6 = 6 0 



T s =T e 



(10) 
(11) 

(12) 
(13) 
(14) 



It is worth mentioning here that the gas phase temperature T y is function of the ambient temperature 7^ and the 
cylinder surface temperature Tp, and the factor a whose value, according to [19], varies between zero for small 
values of Reynolds number and 1 for large values of Reynolds number. 
Thus 



T Y = aT r + (1 - a) 



(15) 



Since the nucleation phenomenon is not well understood, we adopted here that each control volume contains 
only one nucleus and hence we can write 

d 0 = Kf rac 

where Kf rac is the size of the control volume along the direction parallel to the surface and has the values as 
shown in Table 1 where AC is given in kg/m 3 . It is important to mention that 8 0 cannot be zero because of 
numerical reasons, and in the present study a small value given by 8 0 = ^ is adopted. 



Table 1 Adopted value s of Kf 



AC < 0.008 


0.008 < AC 
< 0.01 


0.01 <AC 
< 0.013 


AC > 0.013 


0.5 


0.65 


0.7 


0.75 



The frost formation model divided the process artificially into two stages, but the real formation process is a 
continuous one and the ice column are transformed into porous mesh of different structures. Hence a set of 
coupling equations are necessary to adequate the results of the first stage of column structure as entries to the 
second stage of frost formation. 

These coupling equations are: 



The volumetric fraction £p t is given by 

2 



The temperature 7} can be evaluated by 
Tt=Ep ti T p +(l-Ep ti )T Y 

The density pf t and the specific heat c p ^are evaluated by 
Pf,i = £p,iPp + (1 - £p,dPy 

c p,i = e p,i c p,p,i + (1 — £ p ,i) c p,y,i 



(16) 
(17) 

(18) 
(19) 
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1.2 The crystals spreading phase 

As mentioned before, in this phase the crystals are spread over the cold surface forming a porous medium. This 
model is subject to some restrictions such as heat and mass transfer are one dimensional; 

constant internal pressure in the porous medium; 
thermal equilibrium among the phases; 
dominant thermal diffusion in relation to convection; and, 
Saturated air in the porous medium. 

Considering a control volume in the porous medium, as shown in Fig. 3 one can write the energy balance in the 
form 



dT . - d ( dT\ 



(20) 



The mass balance for a control volume 



Caminhoarbitrario 
De difusao 



Fase solida B 



Fase y 
(vapor+ ar 




or 

Amostra volumetrica 



* * s * . 

.xxxxxxn\>xxxx>.\.xnS\n\\ 

. * h\ xxxxxx\XXXXXXX\xxx\ 
/ 'S^J / 

vxvvvvxv* ' \ • •*XSXX\\\X . 

'.K m' n tZ'^'SS ***** 



,\\\\\XXX '"^nff ' LJa.H \ \ \ \ % \ * 1 V- 

\X\\xx\\\\>\\\x\XX\xxx f 

x Nr.... xxx\\\\\\\\\xx\x \ I 

✓ / f~rOSl ***+*%***%*%*%*%*%*%*%*%*%*%*%* In 



6,(1) 



Fig. 3 Control volume representation in the porous medium after [19]. 
Phase y 

~ dr \ u *f ar ) 



at 



Phase p 

-J- = rn = W -W sc 



(21) 



(22) 



Equations (20) to (22) represent the second stage model. By using adequate initial and boundary conditions 
these equations can be solved. 

It is important to mention that for the solution the equations (20) to (22) the coupling equations (16) to (19) are 
used to adapt the results of the first stage to be used as entry values for the second stage. 
The boundary conditions are: 
On the surface of the cylinder r = 0, 
T s =T e 

or 

The conditions at the moving frost interface at 
r = 8 

h m (W» - W 8 ) = D eftS ^ + p ftS % (23) 
h(T x -T s ) = k ef d £-p f>s A d £ (24) 

^r = ° 

or 

The relevant thermal properties of the frost and air are obtained from the literature as below: 
• The thermal conductivity k Q [24] 

k g = y[W/mK], where T in K. (25) 
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• The specific heat c pg [24] 

c Pt9 = 2116.56 + 7.284510'- 3 7 [kJ/kgK], where T in K. 

• Specific volume v g [25] 

v g = 1.0907xl0~ 3 + 1.4635xl0- 3 7[m 3 /kg], where Tin K. 

• The density p g [16] 

— = 0.0010907 + 1.4635xl0- 7 7[kg/m 3 ], where T in K. 

Pg 

The properties of dry air are obtained from [25]. 



(26) 



(27) 



(28) 



• The thermal conductivity k ar 
k nr = 0.001968 + 8.15xlO- 5 T 



[W/mK], where T in K. 



• The specific heat c par 

Cp ar = 1.004 [kJ/kgK], where T in K 

• The density p ar 

344 9 

Par — — ~ [kg/m 3 ], where T in K. 



The effective diffusion coefficient is evaluated using [19] 
D ef = e y D(l + F)[m 2 /s] 

where D is the binary diffusion coefficient evaluated from 
D = 2.19(0' 10- 5 [m 2 /s] 

• The effective diffusion coefficient on the surface is 
D eff =e Y D(l + F s ) [m 2 /s] 



~i — ■ — i — 1 — i — 1 — i — r 



i — 1 — i — 1 — i — 1 — r 



■ F s = -15.321 + 0.057T 

• F 6 =1 347.83965 + -1 0.72443.T + 0.021 33.T 2 




252 254 256 258 260 262 264 266 268 270 

Temperature [K] 

Fig. 4 Factors F and F s as a function of temperature, [19]. 
where F and Fs are given by [19] and shown in Fig. 4. 



(29) 
.(30) 



(31) 
(32) 

(33) 



(34) 



• The saturation humidity of the air is 
W = 0.6218 ?vav _' sat = 

P total Pvap ,sat 



(35) 



• The saturation pressure at any position in the frost layer is evaluated from [17], 
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ln(P rap ,sat )=f+C 2 + C 3 +C 4 T 2 + C 5 T 3 + C 6 7 4 + C 7 lnRT) (36) 

where 

CI =-5.674 535 9 10 3 
C2 = 6.392 524 7 
C3 = -9.677 843 0 KT 3 
C4 = 6.221 570 1 10" 7 
C5 = 2.074 782 5 10" 9 
C6 = -9.484 024 0 10- 13 
C7 = 4.1635019 



Equation (36) is valid for the interval -100°C<T<0°C where T is in K and P vap sat is in Pa 



• The total specific volume is, 



v = 1.611168. 10" 2 ^(l + Pj f-) [m 3 /kg] (36a) 

• The value of the sublimation latent heat, [12] 

A = h vapJ - h ice>T = -0.195386. T[k] + 2837.77 [kJ/kg] (37) 

The effective thermal conductivity is evaluated by using the model due to [24], where the frost is treated as a 
porous medium of volumetric fraction £ given by 

e = P -^ " (38) 

Pg-Pa 

The effective thermal conductivity is evaluated by 

k ef = \{0£c - Dfci + Ox c - l)k u + V[(3£ c - + (3x c - l)/cj 2 + 8 k t k u ] (39) 
Where 

k u = (1 — s)k h + sk c (Upper limit) 
&! = (! — s)k p + sk s (Lower limit) 



Let us define 
h 

^e.ar 

a= ir~ 

Tck 1? 

k = k 4- 

1 

7i- 



where 

(a) Contribution due to air bubbles (/c fe ) : 



1 + e 



\2 + aJ 



(b) Contribution due to the frost cylinders (fc c ): 
k c = (1 — e)k g + £k ear 

(c) Contribution due to the frost spheres (fc s ): 

3 + e(q - 1) 
3 - £ {—) 

(d) Contribution due to the flat plate frost (k p ) : 
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"-g "-e.ar 

(1 - s)k e ar + sk g 



The value of s c is given by the correlation obtained from experimental data: 

s c = i3.6fe - eoce - s.y f i - §^ + + 9r (£ ":; )2 l 

C W l^V U I 3e 3 - £l e 2 -ei 2(£ 3 -£i)(£ 2 -£i)J 

Where 

^ = °- 1726 (^16) 
£ 2 = 0.751 

~ (273.16) 



71 



£ 3 = 0.3sin<( — 



(273!l6) 



If s < s t , then e c = 0. The value of x c is defined as: x c = 1 — £ c 



The equations of the first stage, the coupling equations and the equations of the second stage together 
with the respective boundary and initial conditions were discretized by using the finite difference 
approximation. The numerical program was written in C++ and the numerical tests were realized to optimize the 
computational grid and make the results independent of the grid size. The value of dr was varied from 10" 1 



to 10" 5 mm and dt from 10s to 10" 4 s. 



The predicted results were found to compare well with the experimental measurements for dr = 10" mm and dt 
= 10~ 2 s, and hence these values were used for all the present calculations. 



III. EXPERIMENTAL MEASUREMENTS 

In order to validate the model and the numerical predictions, an experimental rig was constructed and 
adequately instrumented, Fig. 5. The experimental rig is composed of a wind tunnel with the test section made of 
acrylic sheet of 10mm thickness where the cold tube was fixed vertically and fed by the secondary fluid circuit. 
The air flow in the tunnel was controlled by varying the frequency of the driving motor, and the working section 
was calibrated such that both its mean velocity and the air flow rate are obtained in terms of the frequency of the 
driving motor. The humidity of the air was changed by using a commercial electric humidifier and the air was 
usually maintained as fully saturated. The circulating secondary fluid flow rate was measured by a calibrated 
orifice plate and its temperature at entry and exit of the cylinder were measured by calibrated thermocouples. 



Condenser 




Test section 

Fig. 5 Experimental rig. 



The dry and wet bulb temperatures were measured before and after the test tube. 

The deposition of frost on the surface of the cylinder was measured by photographing the cylinder by 
means of a high speed digital camera fixed normal to the cylinder axis. A precision measuring scale was fixed 
close to the cylinder and was always photographed together with the cylinder. 

The uncertainty in the measurements and error propagation analysis were performed for all the experiments 
realized and the final values are listed in Table 2. 
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Table 2 Results of the uncertainty analysis 



Quantity 


Uncertainty 


Time 


S r = ±0.015 


Frost thickness 


e = ± 0.5mm 


Temperature 


5 T = ±0.5°C 


Humid air velocity 




Cylinder diameter 


S D€xt = ±0-0020ww 


Exit section of the wind tunnel 


±0.002 lww 
S A =±0.91ww : 


Air volume flow rate 


5 0 =±0.00088 mVs 


Air mass flow rate 


S m =±O.012\63kg/s 


Secondary fluid 





At the beginning of the frost deposition process the temperature gradient between the surface of the 
cylinder and the air stream is big and consequently water droplets present in the wet air stream collide against 
the cold surface of the cylinder forming ice crystals. As the process continues more ice crystals are deposited 
over the cylinder surface forming a porous matrix. As the frost layer gets thicker it's the thermal resistance 
increases and the temperature difference decreases provoking the fusion of the solidified ice crystals and 
changing it to the liquid phase. When this occurs, the thermal resistance decreases and the liquid layer freezes 
and more ice particles are formed on the surface of the cylinder. This process is repeated many times and the 
final result is the formation of a layer of very dense frost or a layer of glacial ice. Figs.6 a,b show images of the 
frost porous layer and the frost glacial layer. 




Fig. 6b Formation of the glacial frost. 
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IV. RESULTS AND DISCUSSION 

As was mentioned before some numerical tests were realized to establish the optimum computational 
grids which make the numerical predictions independent of the grid size. The predicted test results were 
compared with the experimental measurements. It was found that dr = 10~ 3 mm and dt = lO ^gave good 
agreements with the experiments. Therefore these values were adopted for all the calculations. 
Fig. 7 shows the variation of frost thickness with time and as function of the surface temperature of the cylinder. 
As can be seen the lower the surface temperature of the cylinder the thicker the frost layer for the same 
conditions and time. 




0 20 40 60 80 100 120 140 160 180 200 

Time [min] 



Fig.7 Variation of the frost layer thickness due to variation of the cylinder surface temperature for 

V ar = 0.03 m 3 /s andra /sc = 0.028^-. 
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Fig. 8 Variation of the frost deposition velocity due to variation of the cylinder surface temperature for V ar = 

0.03 m 3 /s and m fsc = 0.028^-. 

Fig. 8. shows the variation of the frost deposition velocity with time and with the surface temperature of 
the cylinder. As can be seen, the lower the surface temperature of the cylinder the higher the frost deposition 
velocity. Also one can observe that the high initial frost deposition velocity due to the small initial thermal 
resistance. As the time goes on more frost is deposited the surface of the cylinder causing the increase of 
thermal resistance and consequently reducing of frost deposition velocity. The deposition of frost continues until 
towards the end of the process where the frost layer is so thick and the thermal resistance is too high so that the 
deposition velocity is very low and the frost deposition nearly stops. 
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4 D 
3 5 



* 0.05 m 3 /s 

a 0.07 m^s 

V 0,09 mrVs 

4 1,20 mrVs 




80 1 00 1 20 140 160 180 200 

Time rminl 

Fig. 9 Variation of the frost thickness due to variation of the flow rate of the humid air for Tf SC = —10 0 C and 



m 



fsc 



= 0.028^. 



Fig. 9 shows the growth of the frost layer thickness with time for various wet air mass flow rates. As can be seen 
the increase of the wet air mass flow rate decreases of the frost layer thickness. This increase of the wet air mass 
flow rate increases the shear action on the layer, tearing off liquid particles and carrying it away with the air 
stream. 



0.03m«/s 
0.05m'/s 
0.07m*/s 
0.09m*/s 
I.OOmVs 
1.20m s /s 





120 140 160 180 200 



Time [min] 

Fig. 10 Variation of the frost deposition velocity due to variation of the flow rate of the humid air for Tf sc 
-10°Candra /sc = 0.028^-. 



The effect of the mass flow rate of humid air on the frost deposition velocity can be seen in Fig. 10. Where one 
can verify that high air flow rate reduces the deposition velocity. This effect can be explained by the shearing 
action of the high velocity, hence reducing the frost thickness and consequently the frost deposition velocity. 




0 20 40 60 80 100 120 140 160 180 200 
Time [min] 

Fig. 11 Variation of the frost thickness due to variation of the mass flow rate of the working fluid for Tf sc 
-10°Candra /sc = 0.028^-. 
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Fig. 11 shows the effect of increasing the mass flow rate of the secondary fluid on the frost thickness. It is found 
that increasing the mass flow rate, leads to increase the internal heat transfer coefficient and hence reduces the 
thermal resistance, increases the temperature gradient and the deposition of frost. 




80 100 120 
Time fminl 

Fig. 12 Variation of the frost deposition velocity due to variation of the mass flow rate of the working fluid for 
T fsc = -13 °C and V ar = 0.028^-. 



Fig. 12 shows that increasing the secondary mass flow rate, increases the internal heat transfer coefficient, 
increases the heat removal rate and hence increases the deposition velocity. 
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Fig. 13 Comparison between the predicted frost thickness and the experimental measurements. 



Fig. 13. shows a comparison between predicted frost thickness and the experimental measurements. As can be 
seen the agreement is relatively good indicating that the model and the predicted results can be used for 
preliminary analysis of the frost deposition problem. 

Fig. 14 shows the variation of the frost deposition velocity with time both numerically and experimentally. As 
can be seen the agreement is good over the whole period. 
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Fig. 14 Comparison between the predicted frost deposition velocity and the experimental measurements. 
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Fig. 15 shows a comparison between the predicted frost thickness and frost deposition velocity and the 
experimental measurements. As can be seen the agreement is good. One can observe some abrupt change in the 
inclination of the numerically predicted thickness. This variation is only numerical due to the effect of the 
coupling model. The same can be said in relation to Fig. 16, which shows the predicted frost deposition velocity 
compared to the experimental measurements. One can observe the good agreement and the presence of the slight 
kink due to the coupling model. 
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Fig. 15 Comparison between numerical and experimental frost layer growth. 
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Fig. 16 Comparison between numerical and experimental frost deposition velocity. 
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Fig. 17 Effects of the humid air flow rate on the frost deposition rate for -10°C and 0.028kg/s. 
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Fig. 18 Effects of the humid air flow rate on the frost deposition rate for -10°C and 0.028kg/s. 

Fig. 17 and 18 show the effect of the humid air mass flow rate on the frost thickness and deposition velocity. It 
is found that the increase of the wet air mass flow rate leads reducing the frost thickness and the deposition 
velocity because of the shearing action on the frost surface. 
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Fig. 19 Effects of temperature of the cold surface on the frost deposition rate for 0.03m 3 /s and 0.028kg/s. 
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Fig. 20 Effects of temperature of the cold surface on the frost deposition rate for 0.03m 3 /s and 0.028kg/s. 

Fig. 19 shows the effect of the surface temperature of the cylinder on the thickness of the frost layer. 
As can be seen lowering the surface temperature of the cylinder increases the temperature gradient and hence 
increases the deposition of frost on the surface of the cylinder. As the time goes on the rate of deposition is 
slowed down as can be seen from the decrease of the inclination of the deposition curve. 
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Fig. 20 shows the effect of the cylinder surface temperature on the frost deposition velocity. As can be seen, 
lowering the surface temperature leads to increase the deposition velocity. 




D.DZD D.D3 5 D.D3D DJJ35 DJtiD DJH5 D.DSD D.DS5 D.HD 
■lass fbw rale of Ihe secondary fluid [tpfe] 



Fig.21 Effects of mass flow rate of the secondary fluid on the frost deposition rate for -13°C and 0.03m 3 /s. 
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Fig. 22 Effects of mass flow rate of the secondary fluid on the frost deposition rate for -13°C and 0.03m 3 /s. 



Figs. 21 and 22 show the effect of increasing the secondary fluid mass flow rate. A can be seen increase the 
secondary fluid mass flow rate increases the internal Reynolds number and hence the internal coefficient of heat 
transfer. This leads to increasing the frost thickness and the frost deposition velocity. 



V. CONCLUSIONS 

Based upon the results and discussion, one can make the following conclusions. 

1 . The frost layer thickness and the frost deposition velocity decrease with the increase of the volumetric flow 
rate of the wet air. 

2. The decrease of the cylinder wall temperature increases the frost layer thickness and increases the frost 
deposition velocity. 

3. The increase of the mass flow rate of the secondary fluid circulating into the cold cylinder increases the 
frost deposited thickness and the frost deposition velocity. 

4. The comparisons between the numerical predictions and the experimental results show relatively good 
agreement and thus validating the model and numerical predictions. This allows using the model for 
preliminary evaluation of frost formation. 
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